Magnetism and the Trimeron Bond

A review of progress in understanding the Verwey transition in magnetite (Fe3O4) over the past decade is presented. This electronic and structural transition at TV ≈ 125 K was reported in 1939 and has since been a contentious issue in magnetism. Long range Fe2+/Fe3+ charge ordering has been confirmed below the transition from crystal structure refinement, and Fe2+ orbital ordering and formation of trimerons through weak bonding of Fe2+ states to two Fe neighbors has been discovered. This model has accounted for many spectroscopic observations such as the 57Fe NMR frequencies. The trimeron lifetime has been measured, and trimeron soft modes have been observed. The origin of the first to second order crossover of Verwey transitions in doped magnetites has been revealed by a nanoparticle study. Electronic and structural fluctuations are found to persist to temperatures far above TV and local structural distortions track the bulk magnetization, disappearing at the 850 K Curie transition. New binary mixed-valent iron oxides discovered at high pressure are found to have electronic transitions and orbital molecule ground states similar to those of magnetite.


■ INTRODUCTION
Goodenough's seminal "Magnetism and the Chemical Bond" introduced important concepts such as orbital-based superexchange rules for explaining magnetism in solids. 1 The magnetic behavior of many transition metal compounds was rationalized using these rules, including the ferrimagnetism of the original magnetic material magnetite (Fe 3 O 4 ), and related ferrite spinels. However, magnetite also undergoes a change in properties at T V ≈ 125 K reported by Verwey in 1939 that reveals further electronic complexity at low temperatures. 2 The full low temperature crystal structure was determined in 2012 and revealed direct magnetically driven Fe−Fe bonding interactions within three-center trimeron units. 3 This paper will review progress on the long-running effort to understand the Verwey transition of magnetite made over the subsequent decade.
At ambient temperature, magnetite adopts the cubic spineltype structure with an inverse formal charge distribution Fe 3+ [Fe 2.5+ ] 2 O 4 over tetrahedral A and octahedral B sites, shown throughout as A[B] 2 O 4 . Ferrimagnetic order occurs below the Curie transition at T C ≈ 850 K as there are twice as many up-spins at the B sites as there are down-spins at the A sites. Each cation site has a Fe 3+ 3d 5 S = 5/2 core spin, and rapid hopping of the one extra down-spin electron for every two B sites results in minority-spin-polarized electronic conductivity so all B-sites are structurally and spectroscopically equivalent. The Verwey transition at T V ≈ 125 K, where magnetite undergoes a structural distortion and becomes electrically insulating, is observed in measurements of heat capacity, conductivity, magnetization, and many other properties. Progress made on understanding this transition during the 20th Century is covered in an extensive review by Walz. 4 Verwey proposed that the 125 K transition is driven by an ordering of Fe 2+ and Fe 3+ ions at the B-sites equivalent to localization of the minority spin extra electrons, 2 a phenomenon now known as charge ordering that has been verified in many other oxides. 5 However, initial simple charge ordered models were incompatible with crystallographic data, and a complex lattice distortion to a monoclinic √2 × √2 × 2 superstructure of the cubic room temperature spinel lattice was later identified. 6,7 The supercell has Cc space group symmetry and contains 56 symmetry-unique atoms (compered to three in the cubic Fd3̅ m high temperature cell). The complexity of this acentric superstructure in addition to practical difficulties arising from microtwinning of Cc domains below the Verwey transition hampered single crystal diffraction studies of the low temperature structure. Several partial structure refinements using powder diffraction data with symmetry constraints, 8−10 or Fe Kedge resonant X-ray diffraction studies, 11−13 reported some evidence for charge order during 2001−2011.
A full refinement of the low temperature Cc superstructure of magnetite against microcrystal synchrotron diffraction data recorded at 90 K was reported by Senn, Wright, and Attfield (hereafter the SWA model). 3 Analysis of the local distortion modes of the BO 6 octahedra revealed complex patterns of Fe 2+ / Fe 3+ charge ordering and Fe 2+ t 2g -orbital ordering evidenced by compressive tetragonal Jahn−Teller distortions, as shown respectively in Figure 1a,b. A later ellipsoidal analysis of local coordinations in the SWA model also revealed the charge and orbital ordering features. 14 However, additional structural displacements leading to anomalous shortening of some B−B distances showed that the extra down-spin electrons are not fully localized as Fe 2+ states but are instead spread over linear threesite units where weak magnetically driven Fe−Fe−Fe bonding results in highly structured three-site polarons termed "trimerons" (Figure 1c). The low temperature structure can thus also be described as a network of corner sharing trimerons ( Figure 1d). It is notable that, out of many theoretical predictions made prior to publication of the SWA model, one computational study did correctly predict the charge and orbital ordering patterns within the Cc superstructure of magnetite and also reported some of the trimeron distortions. 15 This paper will review progress on magnetite and the Verwey transition over the past decade, showing how experimental and computational results have been used to test and build upon the charge, orbital, and trimeron orderings and other features of the SWA model and describing new iron oxides that have been discovered to have trimeron-based and related ground states. The review is organized into sections that cover insights into (A) the low and (B) the high temperature states of magnetite (below and above T V ); (C) results for off-stoichiometric and cationdoped magnetites; and (D) discoveries of other iron oxides with trimeron-based and related ground states.

■ RESULTS
A. Low Temperature Magnetite (below T V ). The experimental reproducibility of the SWA model for the Cc superstructure was verified by a subsequent study in which 22 high-accuracy structure refinements using synchrotron X-ray data from three different 10−40 μm grains of magnetite were performed at temperatures from 20 to 124 K. 16 Analysis of the coordinates showed little variation across the models except for small thermal changes at temperatures just below T V .
The low temperature Cc crystal structure is complex and difficult to visualize, and so it is useful to represent the 168 independent shifts in (x, y, z) atomic coordinates as 168 equivalent frozen phonon amplitudes. Only one O atom mode is present in the high temperature cubic structure as a static distortion, and the rest all freeze at the Verwey transition. A total of 80 modes are required for the closest centric description (preserving inversion symmetry) in space group C2/c, and an additional 88 are needed for the full acentric Cc description. The 168 modes belong to four classes; Γ, Δ, X, and W point distortions. The magnitudes of all 168 modes in the SWA model have been analyzed, 3,17 and their thermal variations were also reported. 16 Differences between the amplitudes of centric and acentric branches of Δ, X, and W modes were all found to contribute to the significant off-center atomic distortions in the Cc magnetite structure that can lead to ferroelectric and multiferroic properties. It would be convenient to be able to describe the Cc magnetite structure in terms of a few frozen phonon modes, but no good approximation is yet apparent, although brief details of an attempt at mode parametrization are reported. 18 Further diffraction evidence for charge order in the Cc phase of magnetite has come from a resonant multiwave X-ray diffraction study. 19 The use of three-wave diffraction intensities corrected for self-absorption effects that may have affected earlier studies, giving clear evidence for charge ordering at the Bsites in agreement with the SWA model.
Electronic DFT band structure calculations of the Cc magnetite structure have been reported using the SWA model positions 20,21 or with relaxed coordinates. 22 These have confirmed the reported charge and orbital orderings and show Figure 1. Charge, orbital, and trimeron orders in the low temperature Cc supercell of magnetite, as deduced from the SWA refinement model. 3 (a) Distribution of Fe 2+ /Fe 3+ charge states (blue/yellow spheres). (b) Compressive tetragonal Jahn−Teller distortions arising from orbital order within a single Fe 2+ chain shown as long/short bonds (gray/blue lines) to oxygen atoms (red spheres). (c) Single trimeron unit consisting of three Fe sites with parallel S = 5/2 spins as shown by the up brown-green arrows. Orbital order at the central Fe 2+ site localizes an antiparallel spin electron (small down arrow) in one of the t 2g orbitals which distorts the local structure through elongation of four Fe−O bonds and shortening of the Fe−Fe distances through weak bonding to two Fe neighbors in the same plane, as indicated by the purple arrows. The down spin electron density is approximated by the green ellipsoid. ( that the extra electrons occupy a narrow minority-spin band just below the Fermi level. Real space plots of the electron density show a buildup of charge between Fe atoms that form trimeron units, consistent with a weak bonding effect. 20 Interplay between the orbital order and spin−orbit coupling was found to account for the reported magnetoelectric effect in the Cc structure. 21 DFT calculations have also been used to investigate how well the SWA model accounts for spectroscopic observations of the low temperature magnetite structure. 57 Fe NMR is particularly important as it is the only noncrystallographic technique to have resolved signals from all 24 unique Fe atoms (at 8 A sites and 16 B sites) within the Cc cell. 23 Hyperfine fields from DFT calculations were used to compute the 57 Fe resonance frequencies, 24,25 and these are in excellent agreement with reported values as shown for the B sites in Figure 2. These calculations also support the trimeron description as ref 25 notes "the hyperfine anisotropy data obtained from the DFT calculations support the trimeron concept as the central Fe 2+like ions of the suggested trimerons exhibit significantly larger anisotropy than the end ions ... in agreement with expectations deduced from the description of the electron distribution in the trimerons".
The 57 Fe Mossbauer spectrum has also been simulated using hyperfine parameters from DFT calculations based on the SWA model. 26 Approximation to four sextets was found to give good agreement with Mossbauer data from a high quality single crystal of magnetite. The four signals are in an 8:8:5:3 ratio to account respectively for 8 A-site Fe 3+ ions, 8 B-site Fe 3+ ions, 5 B-site Fe 2+ ions where the extra electron occupies d xz or d yz orbitals, and 3 Bsite Fe 2+ ions with the extra electron in the d xy orbital. The latter group is distinct as the d xy orbitals lie perpendicular to the magnetization in the Cc structure giving rise to lower effective magnetic fields and larger electric field gradients, and they also have distinctive NMR frequencies. 24,25 Further Mossbauer and resonant X-ray experiments have suggested that trimeron direction changes around a fixed central Fe 2+ ion when the easy magnetization axis of the Cc phase is switching by an applied magnetic field. 27 Recent phonon calculations for the Cc structure have shown good agreement with inelastic neutron, X-ray, and nuclear scattering data, revealing strong trimeron− phonon coupling, especially for trimerons oriented parallel to the axes of the monoclinic Cc cell. 28 Dynamics of the low temperature phase of magnetite have been explored using coherent and other light sources. The lifetime of individual trimerons was measured in a pump−probe experiment where the effects of femtosecond laser excitation were followed by soft X-ray diffraction. 29 This found that metallization of the low temperature state of magnetite proceeds in two steps. Initial trimeron destruction takes place in 300 fs, with phase segregation into metallic and insulating regions following on an ∼1500 fs time scale. A full study of the photoinduced phase segregation through optical conductivity measurements was subsequently reported. 30 Soft electronic modes of the trimeron order were recently revealed by low temperature optical pump−terahertz probe experiments. 31 These modes show critical softening and so are associated with the Verwey transition, and they most likely correspond to the sliding of trimerons along their long axes.
The Verwey transition is suppressed at a pressure of 8 GPa as confirmed by changes in elastic constants observed in a high pressure study. 32 Pronounced elastic anisotropies in acoustic waves along the cubic-[110] direction were attributed to the presence of the long Fe−Fe−Fe trimeron axis parallel to this direction. The large shape strain at the Verwey transition makes the low temperature phase sensitive to nonhydrostatic stresses, and twin populations are altered. 33 Uniaxial stresses are found to increase T V initially as twin orientations with higher T V 's become favored. Twinning of the Cc structure is eliminated in small particles, and a study of magnetite nanocrystal size effects showed that the Verwey transition is decreased slightly to T V ≈ 120 K at a 20 nm particle size and is fully suppressed in particles below 6 nm. 34 This demonstrates that the minimum coherence distance for the bulk long-range electronic order is around the length of 10 trimerons.
B. High Temperature Magnetite (above T V ). Above the Verwey transition at T V ≈ 125 K, magnetite has the cubic spineltype structure in space group Fd3̅ m with formal charge distribution Fe 3+ [Fe 2.5+ ] 2 O 4 at ambient temperature. A high temperature powder neutron diffraction study revealed changes in the thermal expansion coefficient and variable oxygen coordinate near 700 K that were attributed to the onset of charge transfer between the tetrahedral A and octahedral B sites. 35 This has been confirmed by recent X-ray spectroscopy measurements which showed that charge transfer from B to A sites, represented by x in the formula Fe 3 + 4 , starts near 330 K and increases up to x = 0.125 at 840 K near T C . 36 Migration of Fe cations from octahedral sites to tetrahedral vacancies was reported at higher temperatures.
A key question has been whether disordered charge, orbital, and trimeron correlations persist in the high-temperature cubic phase. Observation of diffuse scattering just above T V shows that local structural correlations are present, and a single crystal X-ray experiment revealed highly structured diffuse scatter that persists to at least 300 K. 37 This has been corroborated by inelastic scattering studies of the lattice vibrations. Raman studies have shown that changes in vibrational modes associated with the Verwey transition occur from T V up to ∼200 K, 38,39 and an inelastic neutron scattering study up to 293 K found discontinuities in transverse acoustic phonons at T V and a decoupling of electronic and phonon dynamics consistent with slow fluctuations of trimerons in the cubic phase. 40 Anomalous Chemistry of Materials pubs.acs.org/cm Review broadening of Δ and X mode phonons up to at least 293 K was reported from an inelastic X-ray scattering (IXS) study. 41 Resonant IXS (RIXS) has been used to explore electronic excitations of the octahedral Fe cations in cubic magnetite, revealing magnetic excitations driven by polaronic distortions that persist to at least 550 K. 42 Other RIXS experiments have shown that that the orbital components of the magnetic moments are ordered noncollinearly at 300 K, consistent with dynamic distortions associated with polaron formation. 43 However, a RIXS-MLD (magnetic linear dichroism) experiment at 170 K revealed that the polarization dependence of the spin− orbital excitations is incompatible with the purely tetragonal Jahn−Teller distortions of the ideal trimeron quasiparticle (Figure 1c) and suggested that trigonal distortions may be more relevant. 44 Analysis of the PDF (pair distribution function) derived from total scattering experiments has been used to evidence local structural distortions within the cubic phase of magnetite. Room temperature X-ray and neutron PDFs of a nanoparticle magnetite sample were found to be not fitted well by the cubic Fd3̅ m structure, and lower symmetry space groups were used to model the local distortions. 45 The SWA model was used to fit the average degree of local distortion over short, medium, and long-range length scales in an X-ray PDF study covering a wide range of temperatures (90−923 K). 46 The resulting plots in Figure 3 show that long-range structural distortions fall sharply to zero just above T V , while medium range distortions persist up to 250−300 K, which matches Raman observations of modes associated with the electronic order. 38,39 However, short-range structural correlations, on the length scale of an individual trimeron, remain present above T V and decrease to zero near the Curie transition at T C ≈ 850 K, following a similar temperature dependence to the reported bulk magnetization. 47 This also matches the thermal transfer of extra electrons (Fe 2+ states) from octahedral to tetrahedral sites seen by X-ray spectroscopy. 36 The weak bonding Fe−Fe interactions in a trimeron require ferromagnetic alignment of the three core S = 5/2 spins so that the extra minority spin electron can be delocalized over the three Fe ions, as shown in Figure 1c. Hence, magnetization is coupled to local Fe displacements to which the X-ray PDF is particularly sensitive. Fe cation displacements due to Fe−Fe bonding interactions emerging below T C were thus identified as the primary driver of the local structural distortions that give rise to the Verwey transition in magnetite.
C. Doped Magnetites. Many cations can be substituted into magnetite to generate the cubic spinel family of ferrites. Comparison of room temperature X-ray and neutron PDFs for MFe 2 O 4 (M = Mn, Fe, Co, and Ni) spinel nanoparticles showed that the M = Mn, Co, and Ni dopants suppressed the local distortions observed for M = Fe magnetite, 45 consistent with loss of the Fe 2+ states associated with local charge, orbital, and trimeron orders.
Magnetite has a small intrinsic range of nonstoichiometry due to iron-deficiency as Fe 3(1−δ) O 4 up to 3δ ≈ 0.035. Studies of nonstoichiometric and lightly cation-doped magnetites showed that T V is suppressed by doping, and a change from sharp first order to broad second order Verwey transitions was reported around hole doping of x = 3δ = 0.012 48,49 as shown in Figure 4.
The lattice distortion associated with formation of the low temperature Cc state is observed in both first and second order regimes, and no change in phonon spectra between the regimes was found in nuclear inelastic scattering experiments. 50 51 Insight into the origin of the first and second order regimes of doped magnetites has recently been provided by a study of slow oxidation of magnetite nanoparticles. 52 This revealed that the Verwey transition is initially suppressed to a minimum value at T V ≈ 80 K, but on further oxidation recovers to a persistent value of T V = 95 K as shown in Figure 5. This variation demonstrates that the Verwey transition is suppressed not only by the doping Figure 3. Thermal variations of local structural displacements due to electronic fluctuations in magnetite measured from below the Verwey transition at T V ≈ 125 K to above the Curie point at T C ≈ 850 K. 46 Displacements are quantified by the Verwey shift parameter f V which is normalized to the average atomic shift at 90 K in the SWA model. f V was fitted to first/second/third unit cell ranges of the X-ray PDF, which describe short/medium/long-range electronic orders. The first unit cell values show that substantial local structural distortions persist up to T C and closely match the reported variation of the bulk magnetization. 47 This demonstrates that the structural and electronic fluctuations  effect from the added oxygen but also by inhomogeneous strains from the concentration gradient developed between the oxygenrich exterior and oxygen-poor interior of the nanoparticles during oxidation, and this was confirmed by quantitative modeling. Observation of the persistent value of T V = 95 K close to the T V ≈ 100 K crossover between first and secondorder Verwey transitions (Figure 4) 47, 48 shows that the crossover corresponds to the intrinsic lower temperature limit of the Verwey transition in homogeneously doped magnetite. Lower T V values down to 70 K in the second-order regime result from additional effects of strain gradients on the transition. Hence, the reported critical doping δ c = 0.0039 at the crossover 48 is identified as the true upper limit for homogeneous oxygen doping of magnetite.
Low temperature structure refinements of several doped magnetites have been carried out using the same microcrystal method as in ref 3. All of these found the same monoclinic Cc lattice distortion as in pure magnetite. Refinements for a Fe 3(1−δ) O 4 material with estimated hole-doping of 3δ = 0.0116, 53 and for a natural mineral sample of composition Fe 2.986 Al 0.007 Si 0.003 Mg 0.002 Mn 0.002 O 4 54 both gave coordinates similar to those of the SWA model, with the same charge, orbital, and trimeron orders apparent, although with some blurring of the local electronic distortions. However, the refinement of a more heavily doped Fe 3−x Zn x O 4 structure with estimated x = 0.0228 found a remarkable suppression of charge, orbital, and trimeron features at one of the eight Fe 2+ sites within the Cc cell. 53 This site is unique in having its trimeron terminated by another Fe 2+ cation, as shown in Figure 1d, and thus was reported as having a lower ionization potential due to electron−electron repulsion. This discovered doping selectivity is remarkable as it corresponds to a "charge order within a charge order" where the rest of the charge, orbital, and trimeron network of magnetite remains robust while one site is preferentially oxidized.
D. Other Iron Oxides. The trimerons observed in the low temperature structure of magnetite are an example of orbital molecules, clusters made up of coupled orbital states on several metal ions within an orbitally ordered (and often also charge ordered) solid. Further examples of orbital molecules are found in other transition metal compounds, e.g., the V 2 dimers formed below the metal−insulator transition in VO 2 , and are reviewed elsewhere. 55 Recent discoveries of trimeron and related dimeron cluster orders in iron oxides are described below.
Fe 3 O 4 was previously the only known stoichiometric, binary, mixed-valent iron oxide, but the past decade has seen an explosion of iron oxide discoveries. The breakthrough occurred when geophysicists exploring possible new iron oxides formed at high pressures and temperatures within Earth's mantle discovered a new composition, Fe 4 O 5 (Figure 6a), that could be recovered to ambient conditions. 56 This has led to discoveries of new binary mixed-valent iron oxides falling into the Fe n O n+1 (n = 4 and 5) and Fe m O m+2 (m = 5 and 7) homologous series. 57 The Fe n O n+1 materials show Verwey-type transitions with charge, orbital, and orbital molecule ordering in their ground states.  Chemistry of Materials pubs.acs.org/cm Review Fe 4 O 5 orders antiferromagnetically below ∼320 K and shows a Verwey-type charge ordering transition at 150 K and a further spin canting transition at 85 K. 58 The Fe 2+ /Fe 3+ charge ordered structure is incommensurate and consists of trimerons and also dimerons. "Dimeron" is used to describe two-center units analogous to trimerons, where one extra electron is shared between two neighboring B site cations with parallel S = 5/2 spins, giving a symmetric Fe 2 5+ dimer. Two further high pressure charge ordered phases, one containing dimerons and trimerons and the second based on dimerons alone, have subsequently been discovered from exploration of the low temperature phase diagram up to 50 GPa. 59 Fe 5 O 6 undergoes a Verwey-type charge ordering transition at 275 K leading to dimeron ordering, with long-range antiferromagnetic order of the dimerons below 100 K, and the dimeron order is reported to be stable to at least 20 GPa. 60 Electronic ordering has also been reported in several ternaries derived from magnetite or the above new iron oxides. High pressure Mossbauer, conductivity, and diffraction studies of the warwickite type oxoborate BFe 2 O 4 (Fe 2 OBO 3 ), which is charge ordered without orbital molecule formation below 280 K at ambient pressure, 61 reported formation of an electron-localized dimeron phase from 16 GPa up to at least 50 GPa. 62,63 CaFe 5 O 7 , the Ca-stabilized version of the as-yet unreported n = 6 member of the Fe n O n+1 family, is reported to show a magnetic transition at 360 K accompanied by a monoclinic lattice distortion, 64 although possible charge and orbital molecule orders are not yet established.
Charge ordering has been reported in MFe 3 O 5 derivatives of Fe 4 O 5 for M = Mn and Ca. Spins in MnFe 3 O 5 order antiferromagnetically below 350 K with further spin transitions at 150 and 60 K. The latter is driven by charge ordering of Fe 2+ and Fe 3+ but without apparent orbital molecule formation. 65,66 CaFe 3 O 5 orders magnetically near 300 K, forming commensurate and incommensurate charge, orbital, and trimeron ordered (CO) phases when stoichiometric, 67 and is notable for displaying electronic phase separation into CO and charge averaged (CA) ground states when slightly doped. 68,69 Direct comparison of local distortions and spin orders from neutronrefined crystal and magnetic structures of coexisting CO and CA phases in a slightly off-stoichiometric Ca 0.96 Fe 3.04 O 5 sample provides a clear demonstration of the conditions for trimeron order (Figure 6b,c). 68 The CO phase has ferromagnetic spin order and Fe 3+ −Fe 2+ −Fe 3+ charge ordering across groups of three edge-sharing FeO 6 octahedra, with Jahn−Teller compression at the central Fe 2+ site and shortening of the Fe−Fe distances. All of these observations are consistent with trimeron formation. However, the CA phase has antiferromagnetic up− down−up spin order across the same sites, with absence of charge and orbital ordering distortions and Fe−Fe shortening, all demonstrating that trimerons are not present.

■ DISCUSSION
The published results described above show that very significant progress in understanding the Verwey transition of magnetite has been made over the past decade. Much of this has been enabled by technique developments. In particular, state-of-the art X-ray synchrotron beamlines have been used for microcrystal determinations of the Cc ground state crystal structure, 3,16,53,54 resonant and multiwave diffraction studies, 19,27 nuclear inelastic scattering, 50 ambient and high pressure powder diffraction studies, 32,45 diffuse and total scattering (PDF) experiments, 37,45,46 X-ray absorption spectroscopy and MCD, 36 and resonant and nonresonant IXS, 41−44 with the majority of these experiments performed at the ESRF. Coherent light sources have also been important for pump−probe studies of lattice dynamics. 29,31 Improved DFT codes have supported many of these investigations, also enabling accurate NMR and Mossbauer spectra to be simulated. [24][25][26]28 Developments in nanoparticle chemistry have enabled influences of particle size and oxygen content on the Verwey transiton to be determined in exquisite detail. 34,52 High pressure discoveries of new mixedvalent binary iron oxides have broadened insights into magnetite. 56,58−60 The SWA model for the Cc crystal structure of magnetite below the Verwey transition has been corroborated by subsequent microcrystal diffraction studies showing only small changes with temperature or doping, although the selective destruction of one trimeron in Fe 2.98 Zn 0.02 O 4 is a notable structural variation. 53 The Fe 2+ /Fe 3+ charge ordering and Fe 2+ t 2g -orbital ordering deduced from the SWA model have been confirmed by many other techniques, notably through spectroscopic assignment of all 24 57 Fe NMR signals and the four classes of Mossbauer resonances. 24−26 The long-running hypothesis that magnetite has a charge ordered ground state, originally proposed by Verwey in 1939, has thus been comprehensively confirmed over the past decade. The origin of the change from first to second order Verwey transitions in doped magnetites has also been revealed by a nanoparticle oxidation study. 52 The trimeron interpretation of the low temperature electronic order has also been supported by subsequent studies; a critical test was again the 57 Fe NMR spectrum where the trimeron model was shown to give a better description than an alternative orbital picture. 25 Measurement of the trimeron lifetime as a distinct step in the photon-induced metallization of magnetite 29 and recent observation of trimeron soft modes 31 add further weight. Comparison of trimeron and nontrimeron ground states in phase separated CaFe 3 O 5 gives direct observation of the conditions for trimeron formation. 68,69 Electronic DFT band structure calculations have confirmed the charge, orbital, and trimeron orderings in magnetite 20−22 and have enabled quantitative interpretation of many spectroscopic results. A recent study (ref 28) concluded that their "results indicate the validity of trimerons (and trimeron−phonon coupling) to explain the physics of magnetite much beyond their original formulation", suggesting that further insights may derive from more sophisticated future theoretical treatments of trimeron quasiparticles.
Understanding of the low temperature state of magnetite has assisted interpretation of the high temperature cubic phase. RIXS has been particularly insightful in showing that charge and orbital fluctuations remain active far above T V . 42 Two studies have revealed similar variations in local distortions 46 and B to A site charge transfer 36 as temperature increases toward T C ≈ 850 K. These paint a consistent picture that long-range magnetic order creates the trimeron bonding distortions that drive charge and orbital ordering, thereby suppressing B to A site charge transfer. Critical fluctuations in the magnetization as temperature increases toward T C thus lead to loss of the local structural distortions and the onset of intersite charge transfer.
However, studies of local structure in the cubic phase of magnetite have not definitively shown that trimeron distortions are responsible for the Fe and other displacements, and one RIXS study suggests that this geometry is not dominant. 44 62,63 show ordering of two-site dimerons or a mix of trimeron and dimeron units, at ambient or high pressure. Hence, the cubic phase of magnetite might contain mixtures of one-site (single Fe 2+ ions), two-site (dimeron), three-site (trimeron), and perhaps other orbital molecule fluctuations that are likely to change their populations with temperature. The nature of magnetite at ambient temperature (above T V ) thus remains a continuing topic for human inquiry, as it has for around 3000 years.

■ CONCLUSIONS
Discoveries over the past decade have led to great progress in understanding of the Verwey transition at T V ≈ 125 K in magnetite. Long range Fe 2+ /Fe 3+ charge ordering below the transition is confirmed from full refinement of the acentric Cc crystal structure, and Fe 2+ orbital ordering and formation of trimerons through weak bonding of Fe 2+ states to two Fe neighbors have been discovered. This model has accounted for many spectroscopic observations such as the 57 Fe NMR frequencies.
The trimeron lifetime has been measured, and trimeron soft modes have been observed. The origin of the first to second order crossover of Verwey transitions in doped magnetites has also been revealed by a nanoparticle oxidation study. Studies of the cubic phase of magnetite have shown that electronic and structural fluctuations persist to temperatures far above T V and local structural distortions track the bulk magnetization, disappearing at the Curie transition at T C ≈ 850 K. However, whether the high-temperature structural fluctuations are trimeron-like remains to be determined. New binary mixed-valent iron oxides discovered at high pressure are found to have similar electronic transitions and orbital molecule ground states, establishing a broader context for the electronic properties of magnetite. The author declares no competing financial interest.

■ ACKNOWLEDGMENTS
The author thanks the coworkers and collaborators shown as coauthors of the research papers reviewed here. Research in Edinburgh was supported by ERC and EPSRC with access to facilities provided by STFC. Prof. J. M. Honig is warmly thanked for the provision of a stoichiometric magnetite sample used in many of our studies.